Chemically-Mediated Quantum Criticality in NbFe2 
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Laves-phase Nbi+cFe2-c is a rare itinerant intermetallic compound exhibiting magnetic quantum 
criticality at Ccr~l-5%Nb excess; its origin, and how alloying mediates it, remains an enigma. 
For NbFe2, we show that an unconventional band critical point (uBCP) above the Fermi level Ef 
explains most observations, and that chemical alloying mediates access to this uBCP by an increase 
in Ep with decreasing electrons (increasing %Nb), counter to rigid-band concepts. We calculate 
that Ep enters the uBCP region for Ccr > 1.5%Nb and by 1.74%Nb there is no Nb site-occupation 
preference between symmetry-distinct Fe sites, i.e., no electron-hopping disorder, making resistivity 
near constant as observed. At larger Nb (Fe) excess, the ferromagnetic Stoner criterion is satisfied. 

PACS numbers: 71.20.Lp,75.10.Lp,75.45.-|-j 



Quantum criticality emerges from collective low-energy 
excitations leading to a second-order phase transition 
at zero temperature (T).G-l In the study of correlated- 
electron materials, e.g., high-T superconductors and 
heavy-Fermion compounds, understanding such critical 
phenomenon remains a principal challenge; and, locating 
any existing quantum critical points (QCPs) is difficult. 
Near these QCPs, quantum fluctuations (rather than 
thermal fluctuations) are observed to give rise to exotic 
effects, e.g., non- Fermi liquid behavior In intermetal- 
lic compounds the situation appears simpler, where, by 
varying non-thermal order parameters, such as applied 
pressure, external magnetic field, or chemical doping, 
various phases near a QCP can be accessed, i.e., ferro- 
magnetic (FM) , antiferromagnetic ( AFM) , and paramag- 
netic (PM) states. Importantly, strong electron-electron 
interactions are not required for an alloy to show corre- 
lated behavior, as now confirmed in Fe-As compounds.!^ 
The Laves phase of Nbi+cFe2_c is one important exam- 
ple and featured in recent reviews of quantum criticality 
in weak magnets.!^ At stoichiometry the susceptibility ex- 
hibits Curie- Weiss behavior down to a spin-density wave 
(SDW) transition at T^dw - 10 K, whereas at --1.5%Nb 
excess a QCP is observed where the SDW collapses and 
non-Fermi-liquid behavior occurs.^ ^ For larger Nb excess 
(c > 0, hole doping) or Fe excess (c <0, electron doping) 
a FM transition is always observed.!^ Indeed, the sen- 
sitivity of the magnetic state to Nb deficiency has long 
been known.'^ 

Although well characterized experimentally, our under- 
standing of the properties of Nbi+cFe2_c is lacking, es- 
pecially how chemical disorder affects the magnetic tran- 
sitions and QCP, and the FM onset at larger dopings. 
Recent theoretical work studied the electronic proper- 
ties of NbFe2,'^ but did not address the critical chemi- 
cal effects. Of course, including properly the effects of 
disorder in these class of systems, specifically at such a 
small doping (c), is a considerable challenge. Attempting 
to address doping, recent studies use the virtual crystal 
approximatiorPSl (VCA) or supercells (ordered array of 
dop ant s),'ii' having severe shortcomings, as we discuss. 



In metals low-energy excitations lie at /near the Fermi 
surface (FS), and, thus, some unique spectral feature 
at/near Ep is t ypical ly required to drive a transition, as 
with FS nestin^i^Iill qj. van Hove band critical point^^ 
(BCP) types of ordering! ^^ ! ^^ Such FS features have been 
rarely identified as the origin of quantum criticality. We 
detail how an unconventional BCP (uBCP), i.e., an ac- 
cidental saddle-point dispersion, above Ep is responsible 
for the QC behavior in metallic Nbi+cFe2-c- We show 
that chemical disorder mediates, via d-state hybridiza- 
tion, access to this uBCP even with fewer electrons for 
increasing %Nb, which rigid-band/ VCA cannot describe, 
effects well known in alloy theory. 

We use an all-electron, Korringa-Kohn-Rostoker 
(Green's function) electronic-structure method and the 
Coherent-Potential Approximation^^ (KKR-CPA) to cal- 
culate the electronic dispersion, density of states, total 
energies and doping site-preferences in Nbi-(_cFe2-c- See 
Ref. [18] for details and examples for (dis)ordered alloys. 
In agreement with experiment, we calculated Ccr "onset" 
at 1.5%Nb, and, at 1.74%Nb, when Ep lies exactly at the 
uBCP, we find that Nb has no site-occupation preference 
between the two symmetry-distinct Fe sites, favoring a 
homogeneous solute distribution and no electron-hopping 
disorder, suggesting a near constant resistivity. We find 
competing FM and AFM states, as observed, that are as- 
sociated with the competing wavevectors from the uBCP. 
We show that the FM Stoner criterion in Nbi+cFe2-c is 
obeyed for larger Nb (or Fe) excess. If the uBCP are 
removed from consideration, none of these results hold. 
We conclude that the QCP occurs via chemical-disorder- 
mediated access to the uBCP above Ep. 

To understand chemical disorder effects, we first need 
to appreciate the structure of NbFe2, which crystal- 
lizes in a CIA hexagonal Laves phase with space group 
PGa/mmc (#194). In terms of crystallography, NbFe2 
EE Nb4Fe^^''^Fe[;^''^ with a 12-atom unit ccU: Fe^^'*) sites 
form two Kagome networks (_L to c-axis) separated by 
Pg(2a) si|;gs jj^ a hexagonal sublattice, and Nb atoms oc- 
cupy the interstices, see Fig. [T]ja,b). The Wykoff posi- 
tions (without inversion) are: Nb at 4f(|,|,a;), Fe'^^"-' 
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FIG. 1. (Color online) (a) Laves unit cell, (b) Fe'®''' sites 
form two Kagome nets, with Nb (gray), Fe*^°' (red) and Fe'^''' 
(blue), with higher (lower) planes shaded, (c) NbFe2 bands 
(Ef at eV). Highlighted box shows an uBCP above Ep. 



at 2a(0,0,0) and Fe^*^^) at 6h(y,2y, f). Calculations 
are performed with measurecP' structural parameters 
a = 4.8401 A, c = 7.8963 A and internal coordinates 
X = 0.0629 and y = 0.1697. Using KKR-CPA we may 
study the effect of (in)homogeneous solute distributions. 
Disorder with antisite Nb on one or both Fe sublattices is 
written as Nb4(Fef_"^\,„)Nb^?,^),)2(Fef_''^\,,,Nb(f,"i,)6, or 
Nb4(Fei_cNbc)8 for the homogeneous case. 

To identify critical electronic features and chemical 
disorder effects, we detail the dispersion e(k; E) and 
density of states (DOS). KKR uses constant-E matrix 
inversion to get e({k};E), rather than constant-k di- 
agonalization to get eigenvalues e(k;{E}). To han- 
dle (dis)ordered cases, we calculate the Bloch spectral 
fmictioiPl Ai3(k;E) on a grid of 32 x 32 x 24 k-points 
to project the dispersion. For ordered cases, AB(k;E) 
yields S(E — e(k;E)), i.e., bands; otherwise it exhibits 
disorder-induced spectral broadening in E and k, related 
to the finite electron scattering length. 

We now show that the c-dependence tied to the ob- 
served quantum criticality arises from the NbFe2 dis- 
persion above Ep. Our bands in Fig. [ijc) agree with 
those from full-potential methods.^ The bands cross- 
ing Ep along F-M arise mainly from Fe^^^^ t2£,-orbitals, 
and lead to saddle-point dispersion slightly (6.6 meV) 
above Ep with an unusual flat dispersive region (uBCP), 
Figs.jljc) and[2];b). These uBCP near Ep are not a result 
of symmetry, but arise from accidental band crossings. 
Experimentally, Crook and CywinskP^l inferred that the 
p^i^h) t2g-orbitals in the Kagome nets plays a critical 
role in the competing magnetic order associated with the 
quantum phase transition (QPT). The Nbi+cFe2-c phase 
diagrarrP shows the QCP onset at ambient pressure at 
Ccr ~ 1.5%Nb (hole-doped) and extends to c ~ 2.0%, 
with FM at larger doping. For a metal, only low-energy 
excitations near Ep can be relevant for such low-T transi- 
tions. By small tuning of a non-thermal order parameter, 
i.e., c, the unusual QPT behavior is observed, and can 
be explained if these uBCPs are accessed. 
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FIG. 2. (Color onhne) (a) AEf(c) shift due to %Nb-excess vs. 
e/a for Nb4(Fei_cNbc)8. Ep lies at uBCPs at Ccr = 1.74% 
but enters spectral tails at 1.65%. Inset (a) shows (un)doped 
AB'"'(k; E) along F-M, and (b) expands around uBCPs (bars 
are spectral widths due to disorder broadening). 



For homogeneous doping, we find a chemical-disorder- 
mediated increase of Ep(c) versus %Nb excess, or 
decreasing electron-per-atom (e/a) ratio. Fig. [2ja). 
Counter to rigid-band concepts, Ep rises to uBCP and 
due to Fe-Nb (bond/antibond) alloying hybridization in 
Kagome nets Fe-bands shift lower but bands from the 
pure Nb-layer remain unaffected. Fig. [2ja) inset. (An 
ordered array of impurities exaggerates the effect, see 
below, showing disorder plays a key role.) At 1.74%Nb 
excess, uBCP lay at Ep, a 6.6 meV shift due to alloy- 
ing and disorder (finite life-time) effects, see Fig. ^h). 
The 6.6 meV (or 77 K) sets the maximum temperature, 
as observed, for these effects to occur on stoichiometry. 
The dispersion and disorder-induced widths along F— M 
for < c < c^P"", Fig. ^h), estimates the QC range. By 
1.65%, Ep enters the spectral tails, giving zero-energy 
excitations into the anomalous dispersion; Ep is maxi- 
mally aligned with the uBCP by 1.74% (giving a Lifshitz- 
type transition, see below), and exits by 2%, where FM 
is observed. We conclude that the QCP occurs from 
alloying/disorder-mediated access to the uBCP inherent 
in NbFe2 dispersion above Ep, detailed more below. 

We have performed supercell calculations to illustrate 
how sensitive AEp(c) is to approximations used to ad- 
dress chemical effects, which ignore disorder. Supercells 
are numerically costly due to large cells needed with de- 
creasing %Nb. We performed calculations at two concen- 
trations with cells constructed by substituting a Nb-atom 
for one Fe-atom in a 2 x 2 x 1 (6.25%Nb excess) and a 
2x2x2 (3.125%Nb excess) supercell. The supercells also 
yield a AEp(c) increase, relative to c = 0%, of 117.1 and 
281.3 meV for 3.125% and 6.25%Nb excess, respectively. 
The CPA shifts are 65.4 and 158.1 meV, respectively, 
which lie on the curve in Fig. [2ja) at smaller e/a. So, 
supercells do not provide the correct occupancy proba- 
bility nor hybridization across the Kagome net, missing 
the key disorder effects. Nonetheless, supercell results 
do reinforce the fact that rigid-band/ VC A concepts are 
invalid, missing the critical alloying and disorder effects. 
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FIG. 3. (Color online) Non-magnetic DOS per formula unit 
relative to Ep for Nbi+cFe2-c at Ccr ~ and 1.74%. Inset: 
DOS at Ccr around Ep (vertical dashed line is Ep at c=0%). 



Figure [3] compares the DOS for Nbi+cFe2-c at Ccr to 
0%, which is similar to full-potential results. '^Ei' Disorder 
broadening for c > 0% is evident. Ep lies near a precipice 
of a DOS depression, which plays a role in forming a 
SDW state at c = 0%. n(EF) is 3.59 states (eV - for- 
mula unit)~^ for NbFe2. Notably, for NbFe2, unlike con- 
ventional BCP,'S1 the saddle-point dispersion at k's asso- 
ciated with the uBCP is (beyond) cubic (efc c>c k^) in one 
direction and quadratic in the orthogonal ky,kz plane, 
yielding an chemically- mediated peak in n(EF) when Ep 
and uBCP are aligned. Fig. [3] in set. From Stoner theory 
with interaction parameter 7p2l a FM instability occurs 
if n^Ep)! > 1. For pure Fe d-electrons, / was reported^^ 
between 0.7 — 0.9 eV. From susceptibility data for NbFe2, 
the Stoner factor [1 — n(EF)/]~^ was estimatecP at ^100. 
Our calculated / for Nbi+cFe2_c is almost constant ver- 
sus c (0.88 < I < 0.9), but n(EF) increases with doping, 
increasing n(Ep)I, e.g., ^(Ef) for 1.74%Nb increases to 
3.90 states-(eV FU)~^. Stoner's criterion is satisfied be- 
yond 2%Nb, as observed, and discussed more below. 

Lifshitz transitions are mediated by FS topology 
changes (e.g., collapse of FS neck or loss of pockets). 
An unconventional Lifshitz transition emerging near a 
marginal QCP was proposed in ZrZn2.'^ Such a zero- 
temperature, pressure-induced transition is also associ- 
ated with access to an uBCP by an increase in Ep. The 
FS topology changes are reflected in a maximum in n(EF) 
at the Lifshitz point, as we found for hole-doping in 
NbFe2 at c^r = 1.74% (Fig. [s] inset). Not only the topol- 
ogy, but the FS volume is strongly dependent on doping, 
and in NbFe2 enhanced with hole doping, and also ob- 
served in a hole-doped Bao.3Ko.7Fe2As2.'^ 

To explore beyond the QCP, we studied doping in the 
Fe- and Nb-rich parts of the phase diagram. The DOS 
(atom- and impurity-projected) for 3.125% Nb- and Fe- 
rich Nbi+cFe2_c is shown in Fig. |4j The Nb-rich (hole- 
doped) n{Fip) increases to 4.66 states-(eV FU)~^. Fe-rich 
(electron-doped) alloys behave opposite to what is ex- 
pected from rigid-band theory; i.e. with electron doping, 
n(EF) rises to 4.75 states-(eV FU)~^, which is mainly 
due to Fe-impurity DOS originating from Nb-layers, see 
Fig. H (lower panel). These values of Nb-rich and Fe-rich 
7i(Ef) satisfy the FM Stoner criterion, as observed.l^ 

We have also studied how the nature of the uBCP (or 
QCP) is affected by inhomogeneous Nb doping on the 



two Fe sublattices. We find that, if only Fe^^'*^ sites are 
doped, the uBCP lies at Ep when cir^^ = 2.34%, for an 
average c™^ = 2.34% x (6/8) = 1.75%. For Fe^'^'^^-only 
case, dispersion is similar to the homogeneous case as in 
Fig. [l|c), because mainly Fe^^''^ t2g-orbitals within the 
Kagome sheets are involved. In contrast, if only Fe^^°^- 
sites are doped, the uBCP lies at Ep when cir°^ — 8.7%, 
i.e., c'^^s = 8.7% X (2/8) = 2.02%. The shift of Ep rela- 
tive to dispersion in this case arises due to hybridization 
of Fe^^^^-states indirectly with Fe^^''-' -states. Dispersion 
with Fe'^'^''^-only doping shows an increased slope of the 
uBCP compared to the other cases, changing the charac- 
ter and c-dependence of the response. 

In Fig. [5] we report Nb site-preference energy 
differences (e('5'')-E(2'^)) versus %Nb-excess in 

^MFe^'S^,m'f,:!MFi'-lie^^^^^ where E(") is 

the total energy when only the Fe^")-sites are doped. 
At cl^P" = 1.74%, where the uBCP lies at Ep, there 
is no energetically favored site occupancy, therefore, 
no electron-hopping disorder effects, explaining the 
observed extremely slow variation of p versus doping.l^ 
With no site preference, a homogeneous solute dis- 
tribution is favored, further supporting our focus on 
this case. These results shed light on the microscopic 
phenomenon occurring at/near the QCP: the alignment 
of Ep with the uBCP at 1.74% provides maximum 
response concomitant with no electron-hopping disorder 
and instability to both AFM and FM fluctuations. 

We now explore the stability of competing magnetic 
phases. Out of several small-cell, magnetic configurations 
in NbFe2, we found a AFM state energetically favored by 
— 10.45 meV/atom over the non-magnetic state, where 
local moments on the Fe'^°'- and Fe^^'^'-sites are aligned 
antiparallel with values of 0.69 /xB and —0.98 pB, respec- 
tively, and the Nb has an induced moment of 0.14 /itB 
aligned with the Fe*^^°^-site . Th e same ground state was 
found in other calculations.'^'^ The phase diagram indi- 
cates a SDW state for NbFe2 with a highly itinerant na- 
ture (as for iron-arsenide^ , which competes with nearby 
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FIG. 4. (Color online) Non-magnetic DOS relative to Ep for 
3.125% Nb-rich (top) and 3.125% Fe-rich (bottom). 
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FIG. 5. Nb site-preference energy difTerence (e'^'^'-E'^"') 
versus %Nb excess in Nb4(Fei_cNbc)8- E'"' is tlie total en- 
ergy when only the Fe^^'-sites are doped. At Ccr = 1.74%, no 
Fe-site is favored, i.e., homogeneous solute distribution. 

FM states. The itineracy of magnetism is clear from the 
spin density on Fe-sites, see Supplement, while carrier 
density peaks at Ccr, Fig. [3] inset. 

To connect to scattering experiments, the uBCP 
e(k) exhibit zero-velocity quasiparticles near Ep at 
symmetry-equivalent, non-special k-points (5o(±l,0) 

and Qo(±5,±^), where Qo ~ 0.25 in units of the F-M 
caliper in Fig. fljc). These Q's remain relevant for off- 
stoichiometric alloys. The susceptibility, i.e., 

x(q)-Ej/(^k)-/(ek+q)][ek+q-ek]-^ , (1) 

can exhibit an enhanced response due to a convol ution 
of (un)occupied states near Ej?, as for FS nesting.E^EH 
Albeit weaker, it can occur from saddle-point topology 
too. From the small caliper of the flat part of the uBCP, 
i.e., 0.055(27r/a), we estimate a |g| of 0.07 close to 
the 0.05 estimated from experiment.'^ 

Thus, from our direct calculations or estimated 



features, PM, FM and SDW (AFM) states compete near 
the QCP, until overwhelmed by a Stoner instability at 
larger |c|, shown above. Such competing magnetic be- 
havior is observed*' near the QCP, where the low-T re- 
sistivity (p(T) ^ T'^) exponent i/ varies between 3/2 and 
5/3, giving unusual non- Fermi liquid behavior. 

Quantifying the anomalous response further requires 
x(q, c, w), which is beyond the present scope. But, the 
change in uj=0 susceptibility, Ax{Q,i^=0), yields the 
change in the DOS at Ep, i.e., An(EF), and, if large, 
a FM instability. Neal et al. have extended the Moriya 
x(q, c — 0,ll!) theory to account for uBCP, which yields 
an anomalous frequency response and competing FM and 
SDW states,^^ agreeing with our direct calculations. 

In summary, we have identified an accidental Fermi- 
surface (non-ideal saddle-point) dispersion as the origin 
for observed behavior associated with the quantum crit- 
icality in Nbi+cFe2-c- We find that Nb (hole) doping 
accesses an unconventional band critical points in NbFe2 
that provide the necessary low-energy excitations for a 
(Lifshitz-type) transitions. This origin explains most of 
the observed doping behavior in this QC intermetallic 
compound, specifically, (i) onset of Ccr %Nb-excess for 
the QCP; (ii) almost constant resistivity versus c, {in) 
competing magnetic states and temperature scale, (iv) 
observed SDW wavevectors, and, finally, (v) the stable 
FM states at large hole- or electron-doping (>2%). To 
explore the nature of the QC transition further, these 
electronic and chemical features can be incorporated into 
a model Hamiltonian, but they would not have been dis- 
covered without the full dispersion and chemical alloying 
effects detailed here. 
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